Whole genome methylation array analysis reveals new aspects in Balkan endemic nephropathy etiology

Background Balkan endemic nephropathy (BEN) represents a chronic progressive interstitial nephritis in striking correlation with uroepithelial tumours of the upper urinary tract. The disease has endemic distribution in the Danube river regions in several Balkan countries. DNA methylation is a primary epigenetic modification that is involved in major processes such as cancer, genomic imprinting, gene silencing, etc. The significance of CpG island methylation status in normal development, cell differentiation and gene expression is widely recognized, although still stays poorly understood. Methods We performed whole genome DNA methylation array analysis on DNA pool samples from peripheral blood from 159 affected individuals and 170 healthy individuals. This technique allowed us to determine the methylation status of 27 627 CpG islands throughout the whole genome in healthy controls and BEN patients. Thus we obtained the methylation profile of BEN patients from Bulgarian and Serbian endemic regions. Results Using specifically developed software we compared the methylation profiles of BEN patients and corresponding controls and revealed the differently methylated regions. We then compared the DMRs between all patient-control pairs to determine common changes in the epigenetic profiles. SEC61G, IL17RA, HDAC11 proved to be differently methylated throughout all patient-control pairs. The CpG islands of all 3 genes were hypomethylated compared to controls. This suggests that dysregulation of these genes involved in immunological response could be a common mechanism in BEN pathogenesis in both endemic regions and in both genders. Conclusion Our data propose a new hypothesis that immunologic dysregulation has a place in BEN etiopathogenesis. Keywords: Epigenetics; Whole genome array analysis; Balkan endemic nephropathy

Quantitative 1-Step DNA Methylation Analysis with Native Genomic DNA as Template

DNA methylation analysis currently requires complex multistep procedures based on bisulfite conversion of unmethylated cytosines or on methylation-sensitive endonucleases. To facilitate DNA methylation analysis, we have developed a quantitative 1-step assay for DNA methylation analysis.

Deep-sequencing identification of the genomic targets of the cytidine deaminase AID and its cofactor RPA in B lymphocytes

The cytidine deaminase AID hypermutates immunoglobulin genes but can also target oncogenes, leading to tumorigenesis. The extent of AID's promiscuity and its predilection for immunoglobulin genes are unknown. We report here that AID interacted broadly with promoter-proximal sequences associated with stalled polymerases and chromatin-activating marks. In contrast, genomic occupancy of replication protein A (RPA), an AID cofactor, was restricted to immunoglobulin genes. The recruitment of RPA to the immunoglobulin loci was facilitated by phosphorylation of AID at Ser38 and Thr140. We propose that stalled polymerases recruit AID, thereby resulting in low frequencies of hypermutation across the B cell genome. Efficient hypermutation and switch recombination required AID phosphorylation and correlated with recruitment of RPA. Our findings provide a rationale for the oncogenic role of AID in B cell malignancy.

No evidence for large differences in genomic methylation between wild and hatchery steelhead (Oncorhynchus mykiss)

When salmonid fish that have been raised in hatcheries spawn in the wild, they often produce fewer surviving adult offspring than wild fish. Recent data from steelhead (Oncorhynchus mykiss) in the Hood River (Oregon, USA) show that even one or two generations of hatchery culture can result in dramatic declines in fitness. Although intense domestication selection could cause such declines, it is worth considering alternative explanations. One possibility is heritable epigenetic changes induced by the hatchery environment. Here, we show, using methylation-sensitive amplified fragment length polymorphism, that hatchery and wild adult steelhead from the Hood River do not appear to differ substantially in overall levels of genomic methylation. Thus, although altered methylation of specific DNA sites or other epigenetic processes could still be important, the hatchery environment does not appear to cause a global hypo- or hypermethylation of the genome or create a large number of sites that are differentially methylated.

Genomic analysis of non-splenic marginal zone lymphomas (MZL) indicates similarities between nodal and extranodal MZL and supports their derivation from memory B-cells

Three distinct categories of marginal zone lymphomas (MZLs) are currently recognized, principally based on their site of occurrence. They are thought to represent unique entities, but the relationship of one subtype with another is poorly understood. We investigated 17 non-splenic MZLs (seven nodal, 10 extranodal) by gene expression profiling to distinguish between subtypes and determine their cell of origin. Our findings suggest biological inter-relatedness of these entities despite occurrence at different locations and associations with possibly different aetiologies. Furthermore, the expression profiles of non-splenic MZL were similar to memory B cells.

Proteins encoded in genomic regions associated with immune-mediated disease physically interact and suggest underlying biology

Genome-wide association studies (GWAS) have defined over 150 genomic regions unequivocally containing variation predisposing to immune-mediated disease. Inferring disease biology from these observations, however, hinges on our ability to discover the molecular processes being perturbed by these risk variants. It has previously been observed that different genes harboring causal mutations for the same Mendelian disease often physically interact. We sought to evaluate the degree to which this is true of genes within strongly associated loci in complex disease. Using sets of loci defined in rheumatoid arthritis (RA) and Crohn's disease (CD) GWAS, we build protein-protein interaction (PPI) networks for genes within associated loci and find abundant physical interactions between protein products of associated genes. We apply multiple permutation approaches to show that these networks are more densely connected than chance expectation. To confirm biological relevance, we show that the components of the networks tend to be expressed in similar tissues relevant to the phenotypes in question, suggesting the network indicates common underlying processes perturbed by risk loci. Furthermore, we show that the RA and CD networks have predictive power by demonstrating that proteins in these networks, not encoded in the confirmed list of disease associated loci, are significantly enriched for association to the phenotypes in question in extended GWAS analysis. Finally, we test our method in 3 non-immune traits to assess its applicability to complex traits in general. We find that genes in loci associated to height and lipid levels assemble into significantly connected networks but did not detect excess connectivity among Type 2 Diabetes (T2D) loci beyond chance. Taken together, our results constitute evidence that, for many of the complex diseases studied here, common genetic associations implicate regions encoding proteins that physically interact in a preferential manner, in line with observations in Mendelian disease.

Population genomic tests of models of adaptive radiation in Lake Victoria Region cichlid fish

Adaptive radiation is usually thought to be associated with speciation, but the evolution of intraspecific polymorphisms without speciation is also possible. The radiation of cichlid fish in Lake Victoria (LV) is perhaps the most impressive example of a recent rapid adaptive radiation, with 600+ very young species. Key questions about its origin remain poorly characterized, such as the importance of speciation versus polymorphism, whether species persist on evolutionary time scales, and if speciation happens more commonly in small isolated or in large connected populations. We used 320 individuals from 105 putative species from Lakes Victoria, Edward, Kivu, Albert, Nabugabo and Saka, in a radiation-wide amplified fragment length polymorphism (AFLP) genome scan to address some of these questions. We demonstrate pervasive signatures of speciation supporting the classical model of adaptive radiation associated with speciation. A positive relationship between the age of lakes and the average genomic differentiation of their species, and a significant fraction of molecular variance explained by above-species level taxonomy suggest the persistence of species on evolutionary time scales, with radiation through sequential speciation rather than a single starburst. Finally the large gene diversity retained from colonization to individual species in every radiation suggests large effective population sizes and makes speciation in small geographical isolates unlikely.

Relationship of porcine IGF2 imprinting status to DNA methylation at the H19 DMD and the IGF2 DMRs 1 and 2

BACKGROUND: Porcine IGF2 and the H19 genes are imprinted. The IGF2 is paternally expressed, while the H19 gene is maternally expressed. Extensive studies in mice established a boundary model indicating that the H19 differentially methylated domain (DMD) controls, upon binding with the CTCF protein, reciprocal imprinting of the IGF2 and the H19 genes. IGF2 transcription is tissue and development specific involving the use of 4 promoters. In the liver of adult Large White boars IGF2 is expressed from both parental alleles, whereas in skeletal muscle and kidney tissues we observed variable relaxation of IGF2 imprinting. We hypothesized that IGF2 expression from both paternal alleles and relaxation of IGF2 imprinting is reflected in differences in DNA methylation patterns at the H19 DMD and IGF2 differentially methylated regions 1 and 2 (DMR1 and DMR2). RESULTS: Bisulfite sequencing analysis did not show any differences in DNA methylation at the three porcine CTCF binding sites in the H19 DMD between liver, muscle and kidney tissues of adult pigs. A DNA methylation analysis using methyl-sensitive restriction endonuclease SacII and 'hot-stop' PCR gave consistent results with those from the bisulfite sequencing analysis. We found that porcine H19 DMD is distinctly differentially methylated, at least for the region formally confirmed by two SNPs, in liver, skeletal muscle and kidney of foetal, newborn and adult pigs, independent of the combined imprinting status of all IGF2 expressed transcripts. DNA methylation at CpG sites in DMR1 of foetal liver was significantly lower than in the adult liver due to the presence of hypomethylated molecules. An allele specific analysis was performed for IGF2 DMR2 using a SNP in the IGF2 3'-UTR. The maternal IGF2 DMR2 of foetal and newborn liver revealed a higher DNA methylation content compared to the respective paternal allele. CONCLUSIONS: Our results indicate that the IGF2 imprinting status is transcript-specific. Biallelic IGF2 expression in adult porcine liver and relaxation of IGF2 imprinting in porcine muscle were a common feature. These results were consistent with the IGF2 promoter P1 usage in adult liver and IGF2 promoter P2, P3 and P4 usages in muscle. The results showed further that bialellic IGF2 expression in liver and relaxation of imprinting in muscle and kidney were not associated with DNA methylation variation at and around at least one CTCF binding site in H19 DMD. The imprinting status in adult liver, muscle and kidney tissues were also not reflected in the methylation patterns of IGF2 DMRs 1 and 2.

Comparison of genomic and proteomic data in recurrent airway obstruction affected horses using Ingenuity Pathway Analysis(R)

BACKGROUND: Recurrent airway obstruction (RAO) is a severe chronic respiratory disease affecting horses worldwide, though mostly in the Northern hemisphere. Environmental as well as genetic factors strongly influence the course and prognosis of the disease. Research has been focused on characterization of immunologic factors contributing to inflammatory responses, on genetic linkage analysis, and, more recently, on proteomic analysis of airway secretions from affected horses. The goal of this study was to investigate the interactions between eight candidate genes previously identified in a genetic linkage study and proteins expressed in bronchoalveolar lavage fluid (BALF) collected from healthy and RAO-affected horses. The analysis was carried out with Ingenuity Pathway Analysis(R) bioinformatics software. RESULTS: The gene with the greatest number of indirect interactions with the set of proteins identified is Interleukin 4 Receptor (IL-4R), whose protein has also been detected in BALF. Interleukin 21 receptor and chemokine (C-C motif) ligand 24 also showed a large number of interactions with the group of detected proteins. Protein products of other genes like that of SOCS5, revealed direct interactions with the IL-4R protein. The interacting proteins NOD2, RPS6KA5 and FOXP3 found in several pathways are reported regulators of the NFkappaB pathway. CONCLUSIONS: The pathways generated with IL-4R highlight possible important intracellular signaling cascades implicating, for instance, NFkappaB. Furthermore, the proposed interaction between SOCS5 and IL-4R could explain how different genes can lead to identical clinical RAO phenotypes, as observed in two Swiss Warmblood half sibling families because these proteins interact upstream of an important cascade where they may act as a functional unit.

The genomic Echinococcus microsatellite EmsB sequences: from a molecular marker to the epidemiological tool

In the field of molecular and epidemiological parasitology, characterization of fast evolving genetic markers appears as an important challenge to consider the diversity and genetic structure of parasites. The study of respective populations can help us to understand their adaptive strategies to survive and perpetuate the species within different host populations, all trying to resist infection. In the past, the relative monomorphic features of Echinococcus multilocularis, the causative agent of alveolar echinococcosis and a severe human parasitic disease, did not stimulate studies dealing with the genetic variability of Echinococcus species or respective populations. A recently developed, characterized and validated original multilocus microsatellite, named EmsB, tandemly repeated in the genome, offered an additional opportunity for this line of investigation. We have compiled in this review new insights brought by this molecular tracker on the transmission activity of Echinococcus among different hosts and at different geographical scales.

Genomic Data Reveal a Complex Making of Humans

In the last few years, two paradigms underlying human evolution have crumbled. Modern humans have not totally replaced previous hominins without any admixture, and the expected signatures of adaptations to new environments are surprisingly lacking at the genomic level. Here we review current evidence about archaic admixture and lack of strong selective sweeps in humans. We underline the need to properly model differential admixture in various populations to correctly reconstruct past demography. We also stress the importance of taking into account the spatial dimension of human evolution, which proceeded by a series of range expansions that could have promoted both the introgression of archaic genes and background selection.

Evolutionary forces shaping genomic islands of population differentiation in humans

Background Levels of differentiation among populations depend both on demographic and selective factors: genetic drift and local adaptation increase population differentiation, which is eroded by gene flow and balancing selection. We describe here the genomic distribution and the properties of genomic regions with unusually high and low levels of population differentiation in humans to assess the influence of selective and neutral processes on human genetic structure. Methods Individual SNPs of the Human Genome Diversity Panel (HGDP) showing significantly high or low levels of population differentiation were detected under a hierarchical-island model (HIM). A Hidden Markov Model allowed us to detect genomic regions or islands of high or low population differentiation. Results Under the HIM, only 1.5% of all SNPs are significant at the 1% level, but their genomic spatial distribution is significantly non-random. We find evidence that local adaptation shaped high-differentiation islands, as they are enriched for non-synonymous SNPs and overlap with previously identified candidate regions for positive selection. Moreover there is a negative relationship between the size of islands and recombination rate, which is stronger for islands overlapping with genes. Gene ontology analysis supports the role of diet as a major selective pressure in those highly differentiated islands. Low-differentiation islands are also enriched for non-synonymous SNPs, and contain an overly high proportion of genes belonging to the 'Oncogenesis' biological process. Conclusions Even though selection seems to be acting in shaping islands of high population differentiation, neutral demographic processes might have promoted the appearance of some genomic islands since i) as much as 20% of islands are in non-genic regions ii) these non-genic islands are on average two times shorter than genic islands, suggesting a more rapid erosion by recombination, and iii) most loci are strongly differentiated between Africans and non-Africans, a result consistent with known human demographic history.